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Abstract 



We propose a microscopic stochastic approach to improve description of 
nuclear dynamics beyond the mean- field approximation at low energies. It 
is shown that, for small amplitude fluctuations, the proposed model gives 
a result for the dispersion of a one-body observable that is identical to the 
result obtained previously through a variational approach. Furthermore, by 
projecting the proposed stochastic mean-field evolution on a collective path, a 
generalized Langevin equation is derived for collective variable, which incorpo- 
rate one-body dissipation and one-body fluctuation mechanism in accordance 
with quantal fluctuation-dissipation relation. 

I. INTRODUCTION 

In the mean-field description of a many-body system, the time-dependent wave function 
is assumed to be a Slater determinant constructed with A number of time-dependent single- 
particle wave functions $j(r*, t). The single-particle wave-functions are determined by the 
Time-Dependent Hartree-Fock equations (TDHF) with proper initial conditions [1-3], 

in^ j (f,t) = h( J >)^ j {f,t) (i) 

where h(p) denotes the self-consistent mean-field Hamiltonian. Rather than the single- 
particle wave functions, in many situations, it is more appropriate to express the mean-field 
approximation in terms of the single-particle density matrix p(f,r f ,t), which is defined as 

p(r, 7, t) = ^(r, t)n^, t) (2) 

3 

where rij denotes occupation factors of the single-particle states. In standard TDHF, there 
are A number of occupied states for which the occupation factors are one, rij — 1, and zero 
for unoccupied states. If the initial state is at a finite temperature T, occupation factors are 
determined by the Fermi-Dirac distribution, rij = l/[exp(ej — p)/T + 1]. The single-particle 
density matrix in the mean-field approximation evolves according to the transport equation, 

ih^p(t) = [h(p),p(t)} (3) 
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The expectation value Q(t) of a one-body observable Q and its dispersion are calculated 
according to, 

Q(*) = E<*i(*)IQI*i(*)>«j ( 4 ) 

and 

= E I < ^WIQI^W > |V(1 - n,-) (5) 

The mean-field approximation includes, so called, the one-body dissipation mechanism and 
therefore it provides a good approximation for the average evolution of the collective motion 
at sufficiently low energies around WMeV per nucleon, at which two-body dissipation and 
fluctuation mechanism do not have an important influence on dynamics. In the mean-field 
approximation, while the single-particle motion is treated in the quantal framework, the 
collective motion is treated nearly in the classical approximation. Therefore, the TDHF 
provides a good description of the average evolution of collective motion, however it severely 
restricts the fluctuations of the collective motion. 

Much effort has been given to improve the mean-field description by incorporating two- 
body dissipation and fluctuation mechanisms [4-6]. The resultant stochastic transport the- 
ory provides a suitable framework for dissipation and fluctuation dynamics of nuclear col- 
lisions at intermediate energies. However, the two-body dissipation and fluctuation mech- 
anism does not play an important role at low energies. Fundamental question is that how 
to improve the mean-field dynamics by incorporating one-body fluctuation mechanism (i.e, 
mean-field fluctuations) associated with the one-body dissipation in a microscopic level? 
In this work, we address this question. We restrict our treatment at low energies at which 
mean-field evolution and the one-body dissipation mechanism provide a good approximation 
and propose a simple stochastic mean-field approach for describing fluctuation dynamics. In 
section 2, we give a brief description of the stochastic approach. In section 3, we illustrate 
that the stochastic approach gives rise to the result, which is identical a previous formula de- 
rive using a variational approach. In section 4, by projecting stochastic mean-field evolution 
on a collective path, we derive a generalized Langevin equation for the collective variable. 
Conclusions are given in section 5. 

II. STOCHASTIC TDHF EQUATION 

In the mean-field framework only source for stochasticity can arise from the initial cor- 
relations. In the standard approach, the TDHF eq. (3) provides a deterministic evolution for 
the single-particle density matrix, starting from a well defined initial state. On the other 
hand, as a result of the correlations, the initial state can not have a well-defined single deter- 
minantal form, but rather it must be a superposition of determinantal wave functions. In the 
mean-field approximation, even though correlations are not propagated, initial correlations 
can be incorporated by considering, not a single initial state but a distribution of initial 
Slater determinants. It is well-known that, such initial correlations can be simulated in a 
stochastic description [7] . In order to develop a stochastic description, we need to determine 
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sufficient number of unoccupied single-particle states in addition to the occupied ones, and 
write the single-particle density matrix in the form, 

p\f, ?, t) = £ $*(r, t; A) p x (t ) $, {f, t; A) (6) 

ij 

Here p^ (t ) are the time-independent elements of density matrix which is determined by the 
initial conditions. The main assumption of the approach is that each matrix element is a 
Gaussian random number specified by a mean value pfj(t ) = SijUj, and a variance 



<H(*o) Spj'Ah) = - rij) + nj(l - m)] (7) 

In these expressions, A represents the event label and 8pfj(t ) represents the fluctuating ele- 
ments of initial density matrix, p^ (t ) = 5ijnj+Sp^(t ). At zero temperature, the occupation 
numbers are one or zero, and at finite temperature rij are determined by the Fermi-Dirac 
distribution. In each event, different from the standard TDHF, the time-dependent single- 
particle wave functions are determined by its own self-consistent mean-field according to, 

ih^$ j (r,t;\)=h(p x )$ j (r,t-,\) (8) 

where h(p x ) denotes the self-consistent mean-field Hamiltonian in the event A. Similar to 
eq. (3), we can express the stochastic mean- field evolution in terms of the single-particle 
density matrix as, 

ih^p\t) = [h(p%p\t)} (9) 

As we illustrate in the following sections, in the stochastic mean-field approach, one-body 
dissipation and fluctuation mechanism is incorporated into dynamics in a manner consistent 
with the quantum fluctuation-dissipation relation, in a microscopic level. 



III. FLUCTUATIONS OF ONE-BODY OBSERVABLES 

In the proposed stochastic approach, employing eq.(9) with random initial conditions, 
we generate an ensemble of different events. Each event evolves according to its own self- 
consistent mean-field potential. In this approach, not only the mean value of an observable 
represented by a one-body operator Q, but we can calculate the probability distribution of 
the observable. The expectation value of an observable Q in an event is determined as, 

Q\t) = J2< Ht\ A)|Q|*i(*; A) > p^to). (10) 

ij 

We note that the influence of the initial fluctuations appears in both the elements of the 
density matrix and also matrix elements of the observable. Even if the magnitude of the 
initial fluctuations is small, the mean-field evolution can enhance the fluctuations, and hence 
events can substantially deviate from one another. The mean value of the observable is 
determined by taking an average over the ensemble generated by numerical simulations, 
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Q(t) = Q x (t). In a similar manner, the variance of the observable is calculated according 
to, 



o* Q (t) = mt) - Q{t)f (ii) 

Here, we consider that the fluctuations are small, p x (t) = p(t)+5p x (t), where Sp x (t) shows 
the small amplitude fluctuations of the density matrix around the average evolution p(t). In 
this case, in expression (10) the correlations between the matrix elements of the observable 
and the initial density matrix can be neglected and the mean value of the observable is 
calculated in the usual manner and given by the standard TDHF result, 

Q(t) = '£<* j (t)\Q\* j (t)>n j (12) 
j 

In this expression, the single-particle wave functions are determined according to the stan- 
dard TDHF equations with the standard mean-field Hamiltonian. 

In order calculate the variance of a one-body observable, we need to determine the fluc- 
tuations of the density matrix 5p x (t). Small fluctuations of the density matrix is determined 
by the time-dependent RPA equation, which is obtained by linearizing the stochastic TDHF 
equation (9) around the average evolution, 



ihpp x (t) = [5h x , p(t)\ + [h{p), 5p x (t)]. (13) 



Here, 5h x (t) = (dh/dp) ■ 5 pit) denotes the fluctuating part of the mean-field Hamiltonian. 
Employing a complete set of stationary single-particle representation, we can express the 
matrix elements of the time-dependent RPA equations as, 

ih?-5p x (t)=J2R ijM (t)5p? k (t) (14) 
ol ki 



where, < &i\5p x (t)\&j >= Sp X j(t) denotes the elements of density matrix and the matrix 
R(t) is given by 



%hW = < > s kj - S a < <S> k \h(p)\<f>j > + 

[< ^ k \dh/dp\^i > Pn3 {t) - p in (t) < Q^ldh/dplQjfy >] (15) 

n 

Using super-space notation, eq.(14) can be written in a compact form as, 

ih^ < ij\6p\t) >=< ij\R(t)\5p x (t) > (16) 

or 

ihj t \5p x (t)>=R(t)\5p x (t)> (17) 

where \Sp(t) > acts like a vector with double indices and R(t) is a matrix in this super-space. 
Because of linear form, formal solution of this equation is given by, 
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\5p x (t) >= exp 



- % - f R(s)ds 
n J t 



\Sp\t ) > 



The fluctuating part of an observable is calculated according to 



5Q\h) =< Q|exp 



h 



to 



R(s)ds 



\5p\t ) > 



(18) 



(19) 



where t\ represent the final time at which the observation is made and t is the initial time. 
It maybe more convenient to define a time-dependent one-body operator B{t) according to 



< B(t)\ =< Q exp 



i f fl 
J R(s)d& 



h 



(20) 



It is easy to show that time evolution of the Heisenberg operator B{t) is determined by the 
dual of the time- dependent RPA according to 



ih^B(t) = [h(p),B(t)]+tr i^j-mip] 



(21) 



The solution is determined by backward evolution this equation with the boundary condition 
B{ti) = Q. In an event A, the expectation value of the observable can be expressed as, 



8Q X {ti) =< B(t )\5p\t ) >=]T< Wo)|^-(*o) > S PjM 



(22) 



Then, the variance of the observable is calculated as, 



*q(*i) = E < Uto)\B(t )\^(t ) >< <S> f (t )\B(t )\Mto) > Sp^o) Sp$ f (to) 

= J2\< Ht )\B(to)\^j(t ) > \ 2 n t (l -nj) (23) 

where the last equality is obtained using the formula (7). This result is identical with the 
formula derived in a previous work employing variational approach [8]. 



IV. PROJECTION ON COLLECTIVE PATH 

In order to illustrate the fact that the stochastic TDHF equation describes dynamics 
of fluctuations in accordance with the one-body dissipation mechanism, we give another 
example in this section. We consider that the collective motion is slow and maybe describe 
by a few relevant collective variables. For example, in induced fission dynamics relevant 
collective variables maybe taken as the relative distance of fragments, mass-asymmetry and 
neck parameter. Here, we consider a single collective variable q(t), and introduce the quasi- 
static single-particle representation, 

h(q)y j (r;q) = e j (q)y j (f;q) (24) 

where h(q) = h[p(q)} denotes the mean-field Hamiltonian, in which the time dependence of 
the local density p(r, r; q) is parameterized in terms of the collective variable in a suitable 
manner. We expand the single-particle density in terms of the quasi-particle representation, 
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p(f,f;q) = E*;tf?WO*i(^?)- (25) 

kl 

Both the elements of density matrix pki{t) and collective variable q(t) are fluctuating quan- 
tities. In this section for clarity of notation, we ignore the event label A on these quantities. 
We determine the matrix elements of density in the lowest order perturbation theory in dy- 
namical coupling < ^> k \d^i/dq > q(t). It is preferable that the wave functions are close to 
diabatic structure. Since in the diabatic representation, dynamical coupling is expected to be 
small, hence, it can be treated in the weak-coupling approximation. Diabatic single-particle 
representation can approximately be constructed by ignoring small symmetry breaking terms 
in the mean- field potential [9,10]. 

In order to determine the temporal evolution of collective variable, we use the total 
energy conservation, 

E = J2 < * k (q)\T\%(q) > p lk (t) + ^ Pji (t) < Mq)Mq)\V\Mq)^M > pM (26) 

Ik ijlk 

In the many-body Hamiltonian, for simplicity we take an effective two-body interaction 
potential energy V. The total energy depends on time implicitly via collective variable q(t) 
and explicitly via matrix elements pi k (t). Energy conservation requires, 

dE .dE dE n 

*=^ + ar = (27) 

In this expression —dE/dq represents a dynamical force acting on the collective variable. 
The force depends on time, and it evolves from an initial diabatic form accompanied with 
deformation of Fermi surface towards an adiabatic limit associated with the adiabatic po- 
tential energy. The second term represents the rate of change of the energy due to explicit 
time dependence, 

^ = E< > ^Pikit) = Ee fc (?) (28) 

Ik k 

Here, p k {t) = Pkkif) represents the occupation factors of the quasi-static single-particle 
states. It is possible to derive a master equation for the occupation factors by substituting 
the expansion (25) into eq.(9) to give 

d i 

T^pki = -jfkl Pkl + E \Phj < ^j\d^i/dq > - < ^ k \d^j/dq > Pjl ] q(t) (29) 
where e k i = e k {q) — ei(q)- This can be converted into an integral equation as, 
Pki(t) = f dtxY, G kl {tM) [Pkjih) < VjldVi/dq > tl - 

Jto j 

< ^ k \d^j/dq > tl pjiit^qih) +G kl {t,t )pki{to) (30) 

where G k i(t,ti) = exp{—f- f* e k i(q)ds] is the mean-field propagator and p k i(to) denotes the 
initial value of the matrix element. Substituting this expression in the right hand side of 
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eq.(29) and keeping only diagonal elements of density matrix, we obtain a master equation 
for the occupation factors, 



d ft 

ftPkit) = J t dh [Gjk(Mi) < Vkld^j/dq > t q(t) < ^jld^k/dq > tl 4(h) + c.c] (p k - Pj ) - 
[G jk (t, t ) < ^ k \d^j/dq > t q(t)p jk (t ) + c.c] . (31) 

j 

This is a stochastic master equation for the occupation factors. Therefore, it determines not 
only the mean value of the occupation factors, but also their distribution functions. The 
first term on the right hand side is the one-body collision term determined by the dynamical 
coupling, while the last term represents its stochastic part. The stochasticity arises from the 
initial correlations pj k (t ), as described in section 2. As a result, the stochastic part of the 
collision term has a Gaussian distribution with zero mean and a variance specified in terms 
of the variance of the initial matrix elements according to eq. (7). 

Using the master equation for the occupation factors, the rate of change of total energy 
due to explicit time dependence can be expressed as, 



dE 
~dt 



f dh ]T G jk (t, h) [e kj < ^ k \d^j/dq > t q(t) < ^j\d^ k /dq > tl 4(h) + c.c] (p k - Pj ) - 
Jt0 jk 

YHG jk (t,t )e kj < VkldVj/dq > t q(t)p jk (t ) + c.c.]. (32) 

jk 



The magnitude of the coupling matrix elements, < ^/ k \d^j/dq >, should decrease as a 
function of energy difference, mainly as a result of the mismatch of overlap of wave functions. 
The energy dependence of the average behavior of square of coupling matrix element may be 
described by a Gaussian or Lorentzian form factor. In weak-coupling limit, that we consider 
here, decay time of the memory kernel is determined by the correlation time of the coupling 
matrix elements defined as r c = h/A, where A is the energy range of the form factor of the 
coupling matrix elements. During time interval smaller than the correlation time, t — ti < r c , 
we neglect the variation of single-particle energies as a function of collective variable and 
approximate the mean-field propagator by G k i(t,ti) = exp[—jr(t — t\)e k i\. Furthermore, 
for slow collective motion it is possible to carry out an expansion in powers of t — t±, as 
it was done in the linear response treatment of ref. [11]. Here, we do not carry out such 
an expansion, but consider the collective motion in a parabolic potential well . In this 
case, the memory effect can be taken approximately into account by incorporating harmonic 
propagation of collective motion during short time intervals according to, 

g(t x ) = qtycosSlfa -t) + ^g(t)smfi(t! - t) (33) 

where Q is the curvature of the harmonic well. Incorporating the result into eq.(32), and 
factoring out q(t) from each term, we deduce a generalized Langevin equation of motion for 
the collective variable [12,13], 

„,.. ldM . 2 dE . _ . 
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where M, 7 and £ denotes the inertia, the friction coefficient and the stochastic force, 
respectively. These quantities are given by, 

^a'EKW^P M2 +W pk (35) 



7 (n) 



E 

jk 



77 



77 



( 6jfe _^)2 +?7 2 ( e . fc + ^)2 + r? 2 



Pfc(l-Pi) ( 36 ) 



and 



£(*) = E G i*(Mo) < *fc|^-|*j > ft*(fo) 



(37) 



where ejfc = £j — £k- A similar treatment can be carried out for a harmonic potential barrier 
[13]. In obtaining the expressions for the inertia and the friction coefficient, we let the upper 
limit of time integration in eq.(32) go to infinity, t— 1 — > 00. We notice that the principle part 
of the integration determines the inertia and real transitions described by the delta-function 
part of the propagator determines the friction coefficient. We note that the inertial force 
term in eq.(34) for the collective variable arises from the q(ti) dependence on the coupling 
matrix element in eq.(32). This dependence is included by expanding < ^j\d^k/dq >t ± 
around q(t) in the lowest order in q(ti) — q(t) and using harmonic propagation for q(ti) 
according to eq.(33). Also, a contribution to friction arises from this term, which vanishes 
at zero frequency limit. Here, we ignore this contribution to the friction. The expression of 
the inertia is different from the adiabatic mass formula, which is known as Cranking formula 
[1]. Due to factor €jk in nominator, unlike the adiabatic mass formula, large contributions 
to inertia near quasi-crossings vanishes, and hence rapid variations of the adiabatic inertia 
in the vicinity of quasi- crossing of levels do not show up in the expression (35). 

As a result of stochastic properties of the initial correlations, stochastic force has a 
Gaussian distribution with zero mean = and a second moment determined by the 
autocorrelation function, 



£(f)£(*i) = E G *(Mi) 

jk 



< > 



1 - pk) + pk(l - pj)] ■ 



(38) 



We introduce the Fourier transform of the autocorrelation function of the stochastic force. 
After a straightforward manipulation, it is easy to show that the autocorrelation function is 
related to frequency spectrum of the friction coefficient according to 



mm = £7 d ^ {t - ti] 7 (c) 



(39) 



where 7(0;) is given by eq.(36) in which the collective frequency Q is replaced by the running 
frequency uj. This result represents the quantum fluctuation-dissipation relation associated 
with the one-body dissipation mechanism and it naturally emerges from the stochastic ap- 
proach presented in this work [14,15]. In the expression of the friction coefficient involves 
a form factor arising from energy dependence of the matrix elements. As discussed above, 
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the form factor can be taken as a Gaussian or Lorentzian specified by the correlation energy 
range A. The correlation energy acts as a cut-off for the integral over the frequency in 
eq.(39). At high temperature limit, T >> A, since coth(hu/2T) = 2T/hu>, the expression 
(39) reduces the classical result, 

£(*)£(ti) ~ 2T 70 6{t - h) (40) 

where 70 = 7(0) is the zero frequency limit of the friction coefficient. In the zero temperature 
limit, on the other hand, coth(hu/2T) « 1, and the autocorrelation function, 

/+°° Hiii 
(*"*!) hu 7 ( W ) (41) 

purely associated with quantum zero point fluctuations. 

V. CONCLUSION 

In this work, we propose a microscopic stochastic approach to improve description of 
nuclear dynamics beyond the mean-field approximation at low energies. In order to demon- 
strate scientific validity of the approach, we present two illustrations. In the first case, we 
show that, for small amplitude fluctuations, the proposed stochastic model gives a result for 
the dispersion of a one-body observable that is identical to previous result obtained through 
a variational approach. In the second case, by projecting the stochastic mean-field evolution 
on a collective path, we derive a generalized Langevin equation for collective variable, which 
incorporate one-body dissipation and one-body fluctuation mechanism in accordance with 
quantal fluctuation-dissipation relation. We believe that such illustrations provide a strong 
support that the proposed model provides a consistent microscopic description for dynam- 
ics of density fluctuations in nuclear collisions at low energies. In fact, such a stochastic 
approach provides a useful tool not only for nuclear dynamics but also other fields, such 
electron dynamics in atomic and condense matter physics. Numerical simulations of the un- 
constraint 3-D TDHF equations have been carried out with powerful computational tools. 
In the proposed stochastic description, we need to calculate a sufficient number of events 
by numerical simulations. The numerical effort for such mean-field calculations is not very 
extensive and can easily be accomplished with the present day computational power. 
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